from astropy.table import QTable

from ctapipe.irf import (
    RadialEventWeighter,
    spectrum_from_name,
    spectrum_from_simulation_config,
)

table = QTable(
    dict(
        true_energy=[1.0, 2.0, 0.5, 0.2] * u.TeV,
        true_fov_offset=[0.1, 1.2, 2.2, 3.2] * u.deg,
    )
)
weighter = RadialEventWeighter(
    source_spectrum=spectrum_from_name("IRFDOC_ELECTRON_SPECTRUM"),
    target_spectrum_name="CRAB_HEGRA",
    fov_offset_max=5.0*u.deg,
    fov_offset_n_bins=5,
  )

print(weighter.fov_offset_bins)
table_with_weights = weighter(table)